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Abstract 

Inspired by the interpretation of two dimensional Yang- Mills theory on a 
cylinder as a random walk on the gauge group, we point out the existence of a 
large TV transition which is the gauge theory analogue of the cutoff transition 
in random walks. The transition occurs in the strong coupling region, with 
the 't Hooft coupling scaling as a log N, at a critical value of a (a = 4 on the 
sphere). The two phases below and above the transition are studied in detail. 
The effective number of degrees of freedom and the free energy are found to 
be proportional to N 2 ~z below the transition and to vanish altogether above 
it. The expectation value of a Wilson loop is calculated to the leading order 
and found to coincide in both phases with the strong coupling value. 
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Introduction 



Yang-Mills theory in two dimensions is at the intersection of many different fields 
of theoretical physics. It is one example of non trivial completely solvable gauge 
theory [1,2], in which both perturbative and non perturbative effects can be studied. 
Its large N expansion has been proved to describe a two dimensional string theory 
[3,4], namely a theory of branched coverings on a two dimensional Riemann surface. 
Non trivial topological sectors in the unitary gauge also seem to be related to matrix 
string states [5]. Its partition function on a torus can be described in terms of a 
gas of free fermions [10, 11, 12, 13], and the kernel on a cylinder by the evolution 
of a system of N free fermions on a circle, namely by the Sutherland model. 

Some new connections between two dimensional gauge theories and statistical 
mechanical systems were pointed out in [6] where two dimensional gauge theories of 
the symmetric group S n in the large n limit were investigated. Gauge theories of S n 
also describe n-coverings of a Riemann surface and hence they are closely related 
to two dimensional Yang-Mills theories; the relation being essentially provided by 
Frobenius formula that relates U(iV) characters (in a representation with n boxes 
in the Young diagram) to the corresponding characters of the S n group. 

It was shown in [6] that the partition function of a gauge theory on a disc 
or a cylinder can be interpreted in terms of random walks on the gauge group, 
whose initial and final positions are the holonomies at the ends of the cylinder 2 
and the number of steps is the area of the surface. Although the focus in [6] was 
on the discrete group S n the argument can be trivially extended, as it is shown 
in the Appendix of the present paper, to continous Lie groups. Similar results 
were independently obtained in [16]. It is well known in random walks theory that 
after a certain number of steps the end point of the walk becomes independent of 
the starting point: the walker has lost any memory of the point he started from. 
The critical number of steps after which that happens can be exactly calculated in a 
number of situations, and the corresponding transition is known as cutoff transition. 
Given the correspondence between random walks on the group and gauge theory on 
a cylinder, one expects to find the cutoff transition also in gauge theories. Indeed it 
was found in [6] that for an S n gauge theory where the holonomy on each elementary 
plaquette is given by a single transposition 3 , a cutoff transition occurs in the large n 
limit when the number of plaquettes (and hence the area) is \n log n, in agreement 
with previous results in random walks [14]. Models with more general Boltzmann 
weight for the plaquettes have a richer structure of phase diagrams [6] . The stringy 
interpretation of the cutoff transition in the S n gauge theory is the following: beyond 
the transition the string world sheet is completely connected in the large n limit, 
while before the transition the world sheet consists of a large connected part and 

2 In a disc the starting point is the identity of the group, and on a sphere both starting and 
ending points are the identity. 

3 In terms of the string interpretation this means that in each plaquette there is a single 
quadratic branch point connecting two of the n sheets of the world sheet. 
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of a small fraction (in fact vanishing in the large n limit) of disconnected parts. 

Another well known correspondence relates random walks with random graphs 
[15], that is graphs obtained by randomly connecting n points with p links. These 
can be put in correspondence with random walks on S n made of p steps, each step 
consisting of a simple transposition. Two types of transitions are known in the 
large n large p double scaling limit in random graphs: a percolation transition at 
a critical value j3 = f3 c when p = f3n and the cutoff transition at a = 1/2 when 
p = an log n. Beyond the transition, namely for p > l/2nlogn all the n points are 
connected whereas before the transition a vanishing fraction of disconnected points 
survive. 

It is rather natural at this point to look for a similar transition in the large N 
limit of U(iV) gauge theories. A large N phase transition on a sphere and on a 
cylinder in two dimensional Yang-Mills theories - the Douglas-Kazakov phase tran- 
sition - has been known for quite some time. However this is not a cutoff transition. 
In a cutoff transition the partition function on a disc for instance becomes inde- 
pendent on the holonomy on the border of the disc, and this is not the case in the 
Douglas-Kazakov transition. Besides the Douglas-Kazakov transition occurs at a 
finite value of the 't Hooft coupling whereas from the previous examples it appears 
that the cutoff transition occurs when the area scales as logiV at large N. From 
this point of view the Douglas-Kazakov transition appears more similar to the per- 
colation transition in random graphs, although a precise correspondence is still to 
be found. 

The existence of the cutoff transition in the large N limit of 2D Yang-Mills 
on the sphere when the 't Hooft coupling scales as logiV was proved by a simple 
argument in [6]. The present paper is devoted to study such transition further on 
both the sphere and the cylinder, in order to characterize its phases and give some 
physical interpretation. The paper is organized as follows: in Section 1 we review 
the large N Douglas-Kazakov transition and its physical interpretation. In Section 
2 we introduce the cutoff transition on the sphere and study the phases above and 
below the transition. Section 3 is devoted to the transition on the disc and on the 
cylinder and Section 4 to the calculation of the expectation value of Wilson loops. 



The partition function of a pure gauge theory on an arbitrary orientable two- 
dimensional manifold Ai of genus G, p boundaries and area A has been known 
for many years [1, 2]: 



1 Large N transition 
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The sum runs over all irreducible representations of the gauge group, A is the gauge 
coupling, Xr{gi) is the character of the holonomy gi in the representation r and 
C2(r) is the quadratic Casimir operator in the representation r. A is related to the 
actual area of M. through A = \ 2 NA. 

We consider G = manifolds with at most 2 boundaries; i.e. spheres, discs 
and cylinders. Moreover, we will confine our analysis to the unitary groups U (N) 
and SU(N). A third order phase transition in the large N limit was discovered 
in the case of a sphere by Douglas and Kazakov [7] at a critical value A = ir 2 of 
the rescaled area A. This transition appears to separate a weak coupling (A < n 2 ) 
from a strong coupling (A > n 2 ) regime. These results were generalized to the case 
of a cylinder in [8, 9] where the phase transition was also interpreted as a result of 
instanton condensation. 

The partition function on the sphere can be written as a sum over the set of 
integers n\ > n.2 > ■ ■ ■ > that label the irreducible representations of SU(iV) 
and U(iV) 4 : 

Z = e"^ 2 - 1 ) IE"* " n,) 2 e-^f^K) 2 , (2) 

ni>n2--->njv i<j 

The existence of a phase transition at the critical value A = n 2 can be easily 
derived from (2) by noticing that the partition function (2) is exactly the same 
as the one of a gaussian hermitian matrix model but with the integral over the 
eigenvalues replaced by the discrete sum over the integers rij. The solution of a 
gaussian hermitian matrix model in the large N limit is given by Wigner's semicircle 
distribution law for the eigenvalues: 



where the continuum variables 

x=jj, A(x) = £. 
and the corresponding density of eigenvalues 

P(A) = | (4) 

have been introduced. In (2) however the would-be eigenvalues n» are distinct 
integers and as a consequence the corresponding density p(A) in the large N limit 
is constrained by: 

|p(A)|<l VA. (5) 

Hence the Wigner semicircle solution is acceptable only in the weak coupling phase, 
namely for A < n 2 where the condition (5) is fulfilled. In fact the maximum of |p(A)| 

4 In the case of U(iV) the extra condition njy > must be imposed. 
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occurs at A = 0, it increases with the area A and becomes equal to 1 at A = n 2 , as 
easily seen from (3). The solution in the strong coupling phase A > n 2 was found 
in [7] and is expressed in terms of elliptic integrals. In this phase a finite fraction 
of the eigenvalues condenses, namely the distribution p(A) is flat and equal to one 
in a symmetric interval around A = as shown in fig. 2. 




Figure 1: Eigenvalue distribution for Figure 2: Eigenvalue distribution for 

A < ir 2 A > 7r 2 (i.e. DK transition) 



1.1 Configuration space 

It is well known [10, 11, 12, 13] that the partition function of two dimensional Yang- 
Mills theories with gauge group U(iV) can be interpreted in terms of a gas of N 
free fermions on a circle described by a Sutherland-Calogero model. In particular, 
if we denote by /C 2 (#,0;^4) the kernel on a cylinder of scaled area A and with 
the U(iV) holonomies at the two ends given by the invariant angles 8{ and fa, it 
was shown that K,2{0, fa A) can be interpreted as the propagator in a time A from 
an initial configuration where the positions of the N fermions on the circle are 
given by Q; L and a final configuration with positions labeled by fa. The partition 
function on the sphere of area A is a particular case where the initial and final 
configurations are just 9i — fa — for all i's, namely the amplitude for a process 
where all fermions start at the origin and come back to the origin after a time 
A. By a modular transformation on the kernel of the cylinder one finds that the 
integers rii labeling the irreducible representations of U(iV) are just the discrete 
momenta of the fermions on the circle. While in the momentum representation the 
Douglas-Kazakov phase transition can be interpreted as fermion condensation, in 
the configuration representation it can be seen in terms of instantons condensation 
[8, 9]. In fact, while going from the initial 9i = configuration to the final fa = 
configuration, a fermion can in principle wind an arbitrary number of times around 
the circle. These winding (instantons) configurations do not contribute in the large 
N limit to the weak coupling phase, as shown by the following simple argument. 
Consider the Wigner distribution (3) of momenta in the weak coupling phase. The 
maximum allowed momentum is n max = -t=, hence the maximum shift in position 
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for a single fermion in the time A is given by 



A9 max = A^= = 2VA (6) 

The existence of winding trajectories requires this shift in position to be at least 
2tc, namely A to be greater of tt 2 . Hence the critical value of A, where the Douglas- 
Kazakov phase transition occurs, marks the point where instantons condense, and 
contribute to the functional integral in the large N limit. 

A more detailed understanding of the Douglas-Kazakov phase transition can be 
achieved by introducing in the large N limit the density p(6, t) of fermions in the 
position 6 at a given time (=area) t. Due to the compact nature of the configuration 
space p{6, t) is defined in the interval — n < 9 < ir with p{— ir, t) = p(ir, t). Matytsin 
proved [17] that if the evolution equation of the fermions is given by the Calogero- 
Sutherland model then the density p(6, t) is governed by the Das-Jevicki equation 
[18] which admits Wigner semicircular distribution of radius r{t) as a solution: 



p(9, t) = ^Vr(ty-9 2 , \9\ < r(t) (7) 

with r{t) satisfying the differential equation + = 0. On a sphere of area 
A the boundary conditions are r(0) = r(A) = and the differential equation has 
the solution: 

r(t) = 2] fMEI (8) 

The solution given by (7) and (8) is valid provided the support of the density 
function p(0,t) is in the interval [— at all t, namely provided r(t) < ir. The 
maximum value for the radius r(t) occurs for t = A/2 and is r max = \J~A~. Hence the 
condition for the validity of the Wigner semicircular solution is A < it 2 . Beyond the 
critical value A = ir 2 the fermions "realize" that the space they live in is compact, 
instantons effects become important and (7) is not an acceptable solution any longer. 
On the cylinder a similar phase transition occurs in general, at a critical value of 
the area that depends on the holonomies at the boundaries 5 . If the distribution of 
the invariant angles of the holonomies at the boundaries is the Wigner semicircular 
distribution, then the critical area can be calculated exactly in the same way as for 
the sphere [8]. For a more general discussion see [9]. 



2 Large A(N) phase transition 

It is apparent from the discussion in the previous section and from the explicit form 
(2) of the partition function that as the area of the sphere increases the distribution 

5 However with particular conditions at the boundary the phase transition may also be absent, 
see for instance [21]. 
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of the "momenta" rii (i.e. of the integers that label the U(iV) representations) 
becomes more and more similar to a double step function, like the one drawn 
with dashed lines in fig. 3. In fact for very large areas the attractive quadratic 
potential tends to dominate over the repulsive force produced by the Vandermonde 
determinant. The double step distribution corresponds to the trivial representation 
of U(iV) in which all characters are identical irrespective of their argument. If this 
distribution dominates the functional integral then the kernel on a disc or on a 
cylinder becomes independent from the holonomies at the boundaries. 

As already mentioned in the introduction and discussed in the Appendix, the 
partition functions on the disc and on the cylinder may be interpreted in terms 
of random walks on the group manifold. From this point of view the very large 
area phase, where the sum over the irreducible representations is dominated by the 
trivial representation corresponds to a walk which is so long that the walker has 
lost any notion of the starting point. The transition where this situation sets in 
is known in random walk theory as "cutoff transition", and the same term will be 
used here. 

The cutoff regime occurs for areas larger than a critical N dependent value 
A C (N) which was found in [6]. The argument is very simple: consider the trivial 
double step representation Rq that minimizes the Casimir term £V n 1 

[N-l N-l\ 
R : {ni,--- ,n N } = < — — — } . (9) 



and compare its contribution to the partition function (2) with the one coming 
from a representation R\ in which n\ has been increased by 1, namely in which 
n\ = + 1. In other words we look for the value of A at which R ceases to be 
dominant. A simple calculation shows that the ratio between the two contributions 
is 

Z (A,N,R ) e# 

Z (A,N,R 1 ) W 1 > 

This ratio is larger than 1 (hence Rq dominates and we are in the cutoff phase) if 
A > 4 log N. 

In order to study this phase transition in more detail, it is convenient to parametrize 
the area A by rescaling it with log N: 

A = a\ogN + f3 (11) 

From the previous argument we expect the cutoff transition to occur at the crit- 
ical value a c = 4, separating two distinct phases . So Yang Mills theory on a 
sphere seems to have four phases altogether: the first two at a = separated by 
the Douglas-Kazakov phase transition at (5 — n 2 , the other two when the area is 
logarithmically rescaled with N that are separated by the cutoff transition. 

The aim of this section is to find the saddle point configuration and the free 
energy in the two phases above and below the cutoff point a c = 4: the cutoff phase 
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a > a c has been discussed above and it is rather trivial, but the phase below a c 
appears as an interesting intermediate phase between the strong coupling phase in 
the Douglas-Kazakov transition and the cutoff phase. Hence we shall concentrate 
on this in the rest of the section. 

P(A) 



•N"~ 



Figure 3: Eigenvalues distribution for < a < 4 



Let us consider again the partition function (2) and the corresponding action 



N 



N 



Hi - 7lj 



■\-—y 

ON ^ 



n. 



'12) 



S = 2 J2 ^ 

i>j=l i=l 

We want to find the extremum of this action in the large N limit, when A is 
parametrized as in eq. (11). Since we expect the saddle point distribution of the 
"momenta" rii to be symmetric with respect to the origin ( n« — > — n« ) we shall 
perform the variation only with respect to symmetric configurations, that is we set 



—M <i < M , with M 



N 



and ri- 



(13) 



Using this symmetry one can restrict the sums to non negative values of % {% — 
1, . . . , M) and write the action as: 



M 



M 



M 



S = 2 log (m - n 3 f + Y, log (m + n 3 f - — ^ n\ (14) 

i>j>l i>j>l i=l 

While in the cutoff phase = % for % — 1, • • • , M below the cutoff transition we 
expect a configuration of the type described in fig. 3, namely: 



rii = i i = 1, • • • , M — I 
riM-i+a = M - I + a + r a a — 1, • • • , I 



(15) 
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where the value of I and the spectrum of the integers r a are to be determined. 
With these notations the action can be written as: 



I M-l 



S-Sn 



+ 



a=l j=l 

2£ (log (1 + 



1 + 



a+j 



1 + 



M - I + a + j. 



a + P + 2(M-l) 



) + log (1 + 



a — f3 



— Y,(2(M-l + a)r a + rl), 



(16) 



where So represents the value of the action in the trivial representation R . We 
shall assume that as M — > oo also I — > oo but at a slower rate than M, namely 
l/M — > 0. We shall also assume that I and r a will be of the same order in the large 
M limit. These assumptions will be justified a posteriori, in the sense that they 
will provide a stable saddle point in the large M limit when the area is scaled like 
logM. They are also very reasonable assumptions: I is of order N in the strong 
coupling phase following the Douglas-Kazakov transition and one expects that with 
the logarithmic rescaling of the area it will shrink further by some power of N. 

The first step in dealing with (16) is to make the dependence from M explicit. 
By using the identity 



M-l 



-)(1 + 



j + a' K M - I + j + a 
we can rewrite the action as 



(2M - 21 + a + r a )\ a\ 
[2M-2l + a)\ (a + r a )\ 



(17) 



S-S = 4^1og 



Q = l 



(2M - 21 + a + r a )\ a\ 
(2M-2/ + a)! (a + r a )\ 



+ 2 ^(log (l + - ; £ ±TL — ) + log (1 + 



a + f3 + 2(M -l)- 



a — (3 



A 

2M 



J2(2(M-l + a)r a + r 2 a ), 



(18) 



a=l 



All ratios of factorials in (18) can be reduced to the form log ^^f^ with iV — > oo, 
C 

C — > oo and — — > 0. By repeated use of Stirling formula one finds, up to terms 
C 

that vanish as — — > 0: 



log — C 



logiV + /(^) 



(19) 
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where 



N k=l 



By using this asymptotic behaviour, and introducing continuum variables in the 
large N limit, namely 

x = j r(x) = y J2 = 1 J dx ■ ( 21 ) 

a 

the action finally takes the form 

f 1 M 

S - S = Al 2 dx[r(x)(log(— ) - log + r(x)) + 1 + log 2) + a;(logx 
Jo ' 

- \og(x + r(x)))] +21 2 [ dx [ dy\og(l + r{x) ~ r{y) ) 

Jo Jo x ~y 

- Al 2 i dxr{x) (22) 

Jo 

where subleading terms (by powers of ) have been neglected. 

Let us now parametrize the area A according to eq. (11) and write the action 

as: 

S - S = Al 2 [F log M - F x log I + F 2 \ (23) 
where F , F 1 and F 2 are of order 1 in the large M and / limit and are given by: 



Fo = (1 - £ dxr(x) 

\ — J dxr(x) 
Jo 

f 

Jo 



F 2 = I dx 

+ 



r(x)(— log (x + r{x)) + 1 — log 2 — ^) + a;(logx — log (x + r(x))) 



['dx [ 1 dy\o g (l + r{x) ' r{y) ) 
Jo Jo x-y 



2 Jo 

In order to find the configuration that maximizes the functional integral in the large 
M limit we take the variation of (23) with respect to both I and r(x). The variation 
with respect to I gives the equation: 

(l-f)logM + (j|-i)=logJ (24) 
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which shows that / grows like M 1_ t. This is consistent with what we expected: 
for a — > it gives I ~ M as in the strong phase beyond the Douglas-Kazakov 
transition, and at the cutoff point a = 4 the power vanishes as expected. 
The variation with respect to r(x) gives on the other hand 

1 1 

d y — ; — T~\ ---\og(x + r(x)) = C (25) 

o x + r(x) - y - r(y) 

with 

C = log I - (f - j) log M + t + log 2 = j| - I + t + log 2 (26) 

If one introduces the new variable (, = x + r(x) and the density function p(£) = ^| 
with support in the interval [0, a] equation (25) becomes 6 : 

drjp^--log^C. (27) 

This is a standard type of equation for the density of eigenvalues in the large 
iV limit of matrix models and can be solved by standard analytic methods (for a 
detailed discussion of this equation see for instance [22]). The resolvent function, 
whose discontinuity across the cut gives the density p(£) 5 is given by: 

H(0 = log£ + C - 2 log (28) 
with the additional condition that for large £ 

m = \ + 0(±) (29) 
The corresponding density is given by 

p(0 = -axccos(i/i). (30) 
7r V a 

This solution obviously describes, through the symmetry (13), both the positive 
and the negative region of rii. 

The asymptotic condition (29) gives the two extra equations 

a = 2 , 

C = log 2. (31) 



6 The lowest extreme of the interval is r(0) which is zero by construction, the upper end a is, 
according with the definition of£, a = 1 + r(l). 
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Eq. (30), together with the first of (31) define r(x) completely, although in an 
implicit way. Hence all integrals involved in the definition of F\ and F2 can be 
calculated. The calculation can actually be done analitically and gives: 

1 F 2 1 3 

* = j. i = - 2 -'i < 32 > 

The second equation could have been derived independently from (26) and (31), so 
it constitutes a non trivial consistency check. We can now write explicitely I and 
the free energy F in terms of the area A(M): 

P i> r-\ - a > r A{M) 

I = e"4M 1_ 4 = Me ~ 

F = S-S = 1 -l 2 = 1 -e^M 2 ^ = ±MV^ (33) 

2 2 2 

Beyond the cutoff transition we have instead F = I = 0, as the dominant eigenvalue 
distribution is given by the double step function sketched in fig. 3. 

The interpretation of these results from the point of view of the free fermion 
description is very clear: beyond the cutoff (a > 4) we are effectively in a zero 
temperature situation where all fermions fill the Fermi sea with no holes. Below 
the cutoff instead (a < 4) some excited fermions and the corresponding holes are 
present in proximity of the surface of the Fermi sea both on the positive and negative 
momentum side. The number of fermions above the sea level is given by I in (33). 
The ratio vanishes like M~~ in the large M limit. This distinguishes this phase 
from the strong coupling phase of the Douglas-Kazakov transition, where such ratio 
remains finite , namely the number of fermions above the Fermi sea level is of order 
M. 

In spite of being described in terms of N free fermions, the free energy is pro- 
portional (with the standard 't Hooft scaling) to N 2 , which reflects the original 
number of degrees of freedom in a unitary N x N matrix model 7 . So it is not 
surprising that the free energy becomes proportional to I 2 in presence of I effective 
fermionic degrees of freedom when the area is rescaled by a log M factor. It is as if 
the effective size of the original matrix had shrunk to I x I. A full understanding of 
the reduction of number of degres of freedom from the point of view of the original 
gauge degrees of freedom is still wanted, although some light on it might be thrown 
by the study of the kernels on the disc and the cylinder in the following sections. 

It is apparent from (33) that the number of effective degrees of freedom M 2_ 2 
is the actual order parameter for the cutoff transition. Incidentally, this unusual 
dependence of the number of degrees of freedom upon a, together with the choice 



7 Thc description in the terms of iV fermions follows the integration over the angular variables 
that reduces the matrix model to an integral over the eigenvalues. The ensueing Vandermonde 
determinant makes the wave function describing the eigenvalues antisymmetric. As a consequence 
the total momentum of M left moving fermions is of order M 2 rather than M, that is of the same 
order as the original number of bosonic degrees of freedom. 
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of a large M limit, makes it rather delicate to classify such a transition according 
to standard terminology. 

Let us finally consider the representations of U(iV) and/or SU(iV) that corre- 
spond to the "momentum" distribution pictured in fig. 3 and given in (30). With 
the group SU(iV) 8 this corresponds to a composite representation in the sense of ref. 
[4], whose Young diagram is shown in fig. 4 where the constituent rpresentations 
are denoted by R and S. The rows in the Young diagram of R and S (that coincide 
in this case) have lengths r a and their total number of boxes |r| is given in the large 
N limit by 

|r| = ^r a ~ I 2 J dxr(x) = h 2 (34) 

If the group is U(iV) an arbitrary number of columns of length N can be added 
or subtracted to the Young diagram of fig. 4, and the composite representation 
can be seen as the direct product of the two constituent representations R and S 
with opposite U(l) charges. In fact the integers labeling the representation S are in 
this case negative, which corresponds to changing signs of the invariant angles 
namely to changing U into W . In the strong couplig regime of the Douglas-Kazakov 
transition (A > n 2 ) the partition function is dominated in the large N limit by a 
composite representation of the same type but with the Young diagram of R and S 
made of rows and columns of lengths of order N, rather than I, and a total number 
of boxes of order N 2 rather than I 2 . 



3 Cutoff transition on the disc and on the cylin- 
der 

In this section we are going to consider the partition function on a disc and on a 
cylinder, with fixed holonomies at the boundaries, in the regime where the area A 
is scaled as in (11). We shall show that the cutoff transition occurs also in this 
case at the same critical value of a, except for some particular holonomies at the 
boundaries. This is in analogy to what happens in the case of the Douglas-Kazakov 
transition which was proved to occur on a cylinder in [8, 9], except possibly for 
some special configurations (see [21]) . 

We shall use the standard expression for the partition function on a cylinder 
with holonomies U and V at the boundaries, which is a particular case of (1): 

Z cyl (U,V) = e-^ N2 -VJ2xr(U)Xr(V)e-^n> (35) 

r 

The partition function on a disc is obtained from (35) by taking for instance U — > 1 
and the partition function on the sphere is recovered by taking the double limit 

8 With SU(iV) the term nf in the action should be replaced by YJ 4 nf — ^-(X)i n i) 2 which 
is invariant under n, — > rii + a. However with suitable choice of a the extra term can be set equal 
to zero. 
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N -3 
N - 2 
N - 1 
N I 



R 



T = SR 



Figure 4: Young Tableau for a Composite Representation. 



U — > 1 and V — > 1. The area A in (35) scales as in (11). The sum over the 
representations r is replaced in the large N limit by the saddle point, namely by a 
representation whose "momentum" distribution is of the type described in fig. 3. 
This corresponds to a composite representation, whose constituent representations 
R and S have rows and columns of order I with l/N — > in the large N limit. The 
first step in evaluating (35) is then to calculate the characters in such limit. 

3.1 Characters in the large N, large / limit 

For simplicity, let us consider first the case where only one constituent representation 
is present, which would be the case if only right (or left) moving fermions were 
present above the Fermi sea. This is given by: 

rii = i for i = 1, . . . , N — 1 

rii = N — I + a + r a for i = N — l + a and a — l,...,l (36) 

The dimension A r of this representation can be written in the form 

A TT \ ( n i ~ n j)\ _ TT oj TT (, r a -r() \yr (N-l + a + r a )\ 

' ~ lj V-j) y^ + ^'i 1 / <*-p)U (N-l + a)\ 

= (N-l + a + r a )\ dr_ E r / n 

| r |lll (N-l + a)\ Irl! [ +U{ N ) \ , 1 ' 

where d r is the dimension of the representation of the symmetric group S\ r \ associ- 
ated to the Young diagram r of rows r a . The order of the symmetric group is the 
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total number of boxes in the Young diagram: |r| = J2 a r a . The dependence on N 
in (37) is explicit: the ratio of factorials in (37) is a polynomial in N of degree |r|. 
If /, a and r a are all of order / the coefficient of 7V' r ' _fe in this polynomial is of order 
l k in the large N, large / limit . Hence we can write 



d r 



\r\ 



^ = R"'1 1+ £*4) fc ] (38) 

k 

where the coefficients c& are smooth in the large /, large N limit. If the limit is 
taken keeping l/N finite, as in the strong coupling phase of the Douglas-Kazakov 
transition, all terms at the r.h.s. of (38) are of the same order and cannot be 
neglegted. On the other hand if the double scaling limit is taken with l/N — > 0, 
as in the previous section, then all terms after the 1 are subleading and can be 
neglected. 

The same argument holds if instead of the dimension of the representation we 
consider a character of U(iV). In fact the celebrated Frobenius formula gives: 

^)^E^n(T) <w 

<reS| r | j 

where Xr( cr ) is a character of the symmetric group S\ r \ in the representation labeled 
by the same Young diagram as Xr(U), Sj are the lengths of the cycles in the cycle 
decomposition of a s j = M) an d k a is the number of cycles in a. By taking 
U — > 1 and comparing (38) and (39) one finds that the contribution of * r j^ when 
a consists of k a cycles is ~ l\ r \~ k ° in the large / limit. So if we take the double 
scaling limit where both I and N go to infinity and the ratio l/N goes to zero, all 
terms in the sum over a in (39) are subleading with respect to the one where a is 
the identity. We obtain: 

Xr{u) = ^.(^wT Nlrl[i + ° m)] = Ar (^wT [i + (4o) 

Notice again that in the strong coupling phase of the Douglas-Kazakov transition, 
where in the large N limit l/N is kept constant, all terms in (39) coming from 
different permutations a are of the same order and cannot be neglected. 



3.2 Partition function on the cylinder and special boundary 
conditions 

As a result of previous analysis we find that in the composite representations that 
are dominant in the region < a < 4 the characters Xr(U) depend only from TrU. 
Hence the partition function on the cylinder will depend only on the trace of the 
holonomies on the boundaries. In fact it is apparent from (40) that going from 
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the sphere (U—l) to the cylinder just amounts to a multiplicative factor \^jrj ■ 
In the case of interest however the composite representation contains both chiral 
and anti-chiral component representations (that is R and S of fig. 4) and not 
just one as in the simplified example discussed above. However in the large N 
limit the two component representations are decoupled [4] and the character of 
the composite representation becomes just the product of the characters of the 
component representations: 

/ TvTI \ 2|t-| 

XcomposUU) = Xr(U) X r(U r ) = A?(| — |) (41) 

By replacing (41) into (35) we find 

Z cyl (U, V) = Zarhen ^ — — \) (42) 

where of course the value of \r\ is the one determined by the saddle point equations 
and the equality holds, in the large N limit, only in the regime where the area A 
is scaled as in (11) with a ^ 0. If we set u = |-^r| and v = |-^-| it is almost 
immediate to see that the multiplicative factor at the r.h.s. of (42) is equivalent to 
replace in the action (23) the constant term j3 in the area A with a $(u,v) given 
by: 

P -> $(u, v) = [3 - log u 2 - log v 2 (43) 

The partition function on the cylinder is then the same as the partition function 
on a sphere whose area is obtained from the area of the cylinder by adding the two 
terms logw 2 and log-u 2 . The latter can be interpreted as the areas of the two discs 
necessary to go, in the given momentum configuration, from U = 1 (resp. V — 1) 
to the holonomy at the boundary with |^| = u (resp |^| = v). The areas of the 
two discs are of order 1, so this correction does not affect the position of the cutoff 
transition that remains on the cylinder at a — 4. 

The discussion above relies on the fact that the leading term in Frobenius for- 
mula (39) comes from the identical permutation. However this is not always true: 
if we take the large N and / limit keeping TrU = 9 then all terms in (39) with a 
containing cycles of length 1 would vanish. Assuming that TrU 2 ^ 0, the term in 
Frobenius formula (39) with the highest power on N would then come from permu- 

\r\ 

tations a all made out of cycles of length 2 and would be of order N 2 instead of 
iV' r L Supposing that both the trace of U and of V vanish the coefficient of logM 
in the first term of (22) would be halved. Correspondingly the critical value of a 
at which the cutoff transition occurs would also be halved and become a c = 2 

Let us make this argument general and more quantitative. Let us assume that 
TrU kl 7^ with TrCP = for j < ki and the same for V with k 2 at the place of ki. 

9 This can happen for instance if the holonomy at a boundary has a symmetry of some sort, for 
instance of the type 6 — » 6 + ir, that is preserved through the limiting process. 
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The leading term in (39) will now be: 



where a consists of ^ cycles of length ki 10 From the discussion following (39) we 

desume that in the large / limit ^4^- ~ /' r ' ( ' 1 ^ By using this asymptotic behaviour 
and eq. (44) we can write the partition function on the cylinder in the large N and 
/ limit as in (23), but with F and F\ replaced by the following expressions: 

F ° = ( w 1 + w 2 -i ) l dxr[x) 

1 1 f 1 

Fx = (— + — ) / dxr{x) (45) 



2/c^ 2/^2 jo 

The dependence of F 2 from r(x) is modified, in a so far unknown way, by next to 
leading terms in MfLU 1 w hile the constant parameter (3 is replaced by 

P^P(u kl ,v k2 ) =p-^-\ogu 2 ki -^-\ogv 2 k2 (46) 

where, with obvious notations, u kl = ^ Tr ^ - and v k2 = ^-jf— -■ Although the 
equation for r(x) cannot be derived, the variation with respect to / gives the scaling 
power of / and the cutoff transition point: 



l- 



l ^ m~ 4 (2ir + 2^ ) (47) 
The cutoff transition occurs then at critical point a c (ki, k 2 ) given by: 

=4(^ + ^) (48) 
which generalizes the result on the sphere. 



4 Wilson loops 

In this section we calculate the expectation value of a Wilson loop in a Yang-Mills 
theory with gauge group U(iV) in the large TV limit. The space-time manifold has 
the topology of a two dimensional sphere whose area scales like log N as in eq. (11). 

10 We assume here for simplicity that we take I — ► oo keeping the total number of boxes in the 
Young diagram multiple of fci and fc 2 . 

11 A lot is known on the characters of permutations with cycles all of the same length, so an 
explicit expression for Fi is probably obtainable. 
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We shall follow th approach of Daul and Kazakov [23] and Boulatov [24] who did 
the calculation for constant 

We may think of the sphere of area A as two discs of areas A\ and A2 (A\ + 
A2 = A) sewed along their common boundary and with holonomy on the boundary 
respectively U and W . The Wilson looop is then given by 



W(A 1 ,A 2 ) = (±=TrU) 



1 

Z 



l - did 2 f dU^TrU X i(U) X 2(U^e-^-^ (49) 
Ri,R 2 J 

where d\ and C\ are the dimension and the Casimir operator in the representa- 
tion Ri, referred to the disc of area Ai, likewise for d 2 and C 2 . The quantity 
J dUTrUxi(U)x2(U*) may be either or 1, namely it is 1 when the Young diagram 
of R 2 is obtained by adding one box to the diagram of Ri and otherwise. That 
is, if Ri is labeled by the integers rt\ > n 2 > ■ ■ ■ > un, R2 is labeled by a set 
on integers where one of the rii is increased by one. Daul and Kazakov used this 
property to get rid of one summation and obtained 

mM,A 2) = i E i E(if n(i+^) 



where the sum over i corresponds to all the possible ways of adding one box to the 
diagram. This is not however the whole result, in fact the original expression is 
symmetric under exchange of A 1 and A 2 , so a term with Ai and A 2 exchanged 12 
must be added to (50) and gives: 



Ri t ■> 

e"^r Cl {e^ + e^ n ^j . (51) 

Moreover, eq. (50), as well as (51), is clearly not symmetric under rii <-> — rij, so we 
can not restrict our considerations to rii > any longer; instead we have to consider 
the whole interval —00 < rii < 00. 

Let us first compute W(Ai, A 2 ) in the frozen phase where the sum over Ri is 
dominated by the trivial representation of dimension d± — 1 labeled by rii = i — ^- 2 ^- 

12 This term originates from the fact that / dUTiUxi(U)X2(U^) is different from zero also if the 
representation conjugate to R\ is obtained from the representation conjugate to R2 by adding a 
box in the Young diagram. 
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Vi and R 2 is the fundamental representation of dimension d 2 — N with C 2 — C\ = N. 
By inserting this into (51) one finds: 

W(A u A 2 ) = (e-^+e-%) (52) 

which is, as expected, a typical strong coupling result. If for instance A 1 » A 2 , 
then we have 

W(A ± , A 2 ) ~ (53) 
Since A = Ay + A 2 = a\og N + P, this situation can occur in two ways: 

1) A 2 = p 2 , Ay = alogN + fr: then W ~ e"^ 

2) Ai = cti logiV + Pi, A 2 = a 2 \og N + p 2 with a x > a 2 : then W ~ iV~~?. 

We proceed now to evaluate A 2 ) in the phase before the cutoff, namely for 

a < 4. The sum over representations in (51) can be replaced by the contribution 
of the dominant representation in the large iV limit, calculated on the sphere in 
section 2. The saddle point is unaffected by the presence of the extra term 



En('^)(> + >) 



rij — n 

which is subleading with respect to the action. After replacing in (51) the sum with 
the saddle point contribution, some simplifications occur and we get 

w(a u 4)4En(i+ ^hr) + «*"*) • (54) 

The representation in (54) is of the type given in (15), and the sum over % describes 
all possible ways of adding a box to the Young diagram. However the replacement 
rii — > rii+1 is impossible in the region — I < i < I as the resulting sequence of integers 
would not be monotonic increasing. So the sum over i in (54) can be replaced by a 
sum over a with 1 < a < I. As a matter of fact we must consider only positive a's, 
as adding a box to the adjoint representation amounts to symmetrize with respect 
to Ay and A 2 , and it has been taken already into account. Hence in (54) we must 
replace the index i with a, rii with M — I + a + r a while the index j goes from 
— M to M, namely it goes over both the condensed and the non condensed regions. 
With these substitutions the expression for the Wilson loop becomes: 



W(A U A 2 ) = £i[l-± (l -^)] 

a=l ^ a 



where 



&=M^)-*(i-^) +{ A 2 ^A l} (55) 
£ a = a + r a (56) 
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It is convenient as usual to use in the large N limit the continuum variables x — j, 
£(x) = an d the density function = The crucial part of the calculation is 

the evaluation of Y^p=i 1°S (l + | -xip ) wn i cn can be done following ref. [23]. By 
expanding the logarithm one finds: 

t '° g (1 + ^ H / dVP(v) ^ ~ £ s ? (^ P({f (57) 

The first term at the r.h.s. comes from the k = 1 term of the log expansion and 
can be evaluated using eq. (27), the other terms can be calculated as in [23] and 
give log S1 ° J p . By inserting these results into (55) and using the explicit form of the 
solution (30) one finally obtains (neglecting O(j^) terms): 

W(A U A 2 )=2 J dt S -^-(e-^ + e~%) = + e"^ (58) 

which is exactly the same result as in the frozen phase. The result is not trivial, 
but it was somehow to be expected. Both phases, before and after the transition, 
are strong coupling phases and the expectation value of the Wilson loop should be 
in both of them obtained, to the leading order, by filling the loop with elementary 
plaquettes in the fundamental representation. The effects of the transition are 
expected to appear only at the next-to- leading order ( ~ j^) which is sensitive to 
the O(jj) degrees of freedom which are not frozen below the cutoff. 
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A Gauge theories as random walks 

The equivalence between random walks on a group G and two dimensional gauge 
theories with gauge group G was pointed out in [6]. This equivalence states that 
the partition function of a gauge theory on a disc of area A and holonomy g on the 
boundary coincides with the probability that a suitably defined random walk on the 
group leads, in a number of steps proportional to A, from the identity in the group 
to the group element g. A similar relation holds for the partition function on the 
cylinder. In ref. [6] the attention was focused on discrete groups, where the number 
of steps can be directly identified with the area measured in suitable units, and 
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in particular on the symmetric group S n . In this appendix we extend, in a rather 
straightforward way, the argument presented there to the case of gauge theories on 
continous Lie groups. Similar conclusions have in the meantime appeared in the 
literature [16]. 

Let us consider a random walk on G with the transition probability defined as 
follows: if the walker is in g p G G after p steps, then his position after p + 1 steps 
is obtained by left multiplying g p by an element g, chosen in G with a probability 
t(g). We assume t(g) to be a class function, whose character expansion can then be 
written as: 

t(g)=J2 A rXr(g)t r - (59) 
r 

where A r is the dimension of the representation r. Suppose that the random 
walk starts from an element g G G and denote by K p (g, g ) the probability for the 
walker to be in g after p steps. Of course K (g,g ) is a delta function, namely: 

K (g,go) = Y, A rXr(g- 1 9o) (60) 

r 

Given K p (g',g ), the probability for the walker to be in g after p + 1 steps is 
given by: 

K p+ i(9,9o) = f dg' t(g{g')- l )K p {g\g Q ) , (61) 

where dg' is the Haar measure. By using (59), (60) and (61) it is easy to show by 
induction that K p (g,g Q ) is a class function of gg^ 1 . In fact let us assume that this 
is true for K p (g',g ), namely that K p (g',g ) admits the character expansion 

K p (g',g ) = J2 A rki p) Xr(9 / 9o 1 ) (62) 

r 

Then by replacing (62) and (59) into (61) and performing the integration over 
g' using the characters fusion rules, one finds that K p+i (g, g ) admits a similar 
expansion with 

A> +1 ) = irk® (63) 

As a result K p (g,g ) is given by: 

K p (g,g ) = J2^rXr(gg 1 )t p r (64) 

r 

Since we are dealing with a random walk on a continous Lie group, we expect 
the walk to be a smooth path. This is obtained by letting the number of steps 
p go to infinity and at the same time the length of each step to zero. Each step 
corresponds then to a very small move on the group manifold, that is t(g) is close 
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to a delta function 5(g), namely according to (59), t r ~ 1. We implement this by 
choosing 

t r = e - £h{r) , with e^O (65) 

If we set g — 1 in (64), 

C (r) 

lim lim pe = A and h(r) = — — (66) 

we reproduce exactly the partition function of a gauge theory on a disc of area 
A, thus establishing the desired connection. A random walk with an arbitrary 
transition function h(r) will be related to a generalized Yang-Mills theory with the 
corresponding potential [19, 20] . 

The partition function on a cylinder may be associated to the probability of 
a random walk from a generic point g (in general different from the identity) to 
g. However in order to obtain the correct answer we need to consider a random 
walk not on the group manifold itself, but rather in the space of orbits, where all 
group elements belonging to the same equivalent class are identified. This amounts 
to replacing in (64) g — > h~ 1 g h and integrating over h. This produces a new 
transition probability K p (g,g Q ) which is separately a class function of both g and 
go and is given by: 

K P (9,9o) = J2^(9)Xr(9o 1 K (67) 

r 

This coincides with the kernel on the cylinder, with the parameters identified ac- 
cording to (66). 



References 

[1] B. E. Rusakov, Mod. Phys. Lett. A 5 (1990) 693. 

[2] A. A. Migdal, Sov. Phys. JETP 42 (1975) 413 [Zh. Eksp. Teor. Fiz. 69 (1975) 
810]. 

[3] D. J. Gross and W. I. Taylor, Nucl. Phys. B 400 (1993) 181 [arXiv:hep- 
th/9301068]. 

[4] D. J. Gross and W. I. Taylor, Nucl. Phys. B 403 (1993) 395 [arXiv:hep- 
th/9303046]. 

[5] M. Billo, M. Caselle, A. D'Adda and P. Provero, Nucl. Phys. B 543 (1999) 141 
[arXiv:hep-th/9809095]. 

[6] A. D'Adda and P. Provero, Commun. Math. Phys. 245 (2004) 1 [arXiv:hep- 
th/01 10243]. 



21 



[7] M. R. Douglas and V. A. Kazakov, Phys. Lett. B 319 (1993) 219 [arXiv:hep- 
th/9305047]. 

[8] M. Caselle, A. D'Adda, L. Magnea and S. Panzeri, arXiv:hep-th/9309107. 

[9] D. J. Gross and A. Matytsin, Nucl. Phys. B 429 (1994) 50 [arXiv:hep- 
th/9404004]. 

[10] M. R. Douglas, arXiv:hep-th/9303159. 

[11] J. A. Minahan and A. P. Polychronakos, Phys. Lett. B 312 (1993) 155 
[arXiv:hep-th/9303153]. 

[12] J. A. Minahan and A. P. Polychronakos, Phys. Lett. B 326 (1994) 288 
[arXiv:hep-th/9309044]. 

[13] M. Caselle, A. D'Adda, L. Magnea and S. Panzeri, Nucl. Phys. B 416, 751 
(1994) [arXiv:hep-th/9304015]. 

[14] P. Diaconis and M. Shahshahani, Z. Wahrsch. Verw. Gebiete 57 (1981) 159. 

[15] I. Pak and V. H. Vu, Discrete Applied Math., 110 (2001) 251. 

[16] S. de Haro and M. Tierz, Phys. Lett. B 601 (2004) 201 [arXiv:hep-th/0406093]. 

[17] A. Matytsin, Nucl. Phys. B 411 (1994) 805 [arXiv:hep-th/9306077]. 

[18] S. R. Das and A. Jevicki, Mod. Phys. Lett. A 5 (1990) 1639. 

[19] M. R. Douglas, K. Li and M. Staudacher, Nucl. Phys. B 420 (1994) 118 
[arXiv:hep-th/9401062]. 

[20] O. Ganor, J. Sonnenschein and S. Yankielowicz, Nucl. Phys. B 434 (1995) 139 
[arXiv:hep-th/9407114]. 

[21] S. Zelditch, Commun. Math. Phys. 245 (2004) 611 [arXiv:hep-th/0305218]. 

[22] K. Kostov, M. Staudacher and T. Wynter, Commun. Math. Phys. 191 (1998) 
283 [arXiv:hep-th/9703189]. 

[23] J. M. Daul and V. A. Kazakov, Phys. Lett. B 335 (1994) 371 [arXiv:hep- 
th/9310165]. 

[24] D. V. Boulatov, Mod. Phys. Lett. A 9 (1994) 365 [arXiv:hep-th/93 10041]. 



22 



